A MATHEMATICAL CLUE TO THE SEPARATION PHENOMENA ON 
THE TWO-DIMENSIONAL NAVIER-STOKES EQUATION 
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Abstract. In general, before separating from a boundary, the flow moves toward reverse 
direction near the boundary against the laminar flow direction. Here in this paper, a clue 
to such reverse flow phenomena (in the mathematical sense) is observed. More precisely, 
the non-stationary two-dimensional Navier-Stokes equation with an initial datum having 
a parallel laminar flow (we define it rigorously in the paper) is considered. We show that 
the direction of the material differentiation is opposite to the initial flow direction and 
effect of the material differentiation (inducing the reverse flow) becomes bigger when the 
curvature of the boundary becomes bigger. We also show that the parallel laminar flow 
cannot be a stationary Navier-Stokes flow near a portion of the boundary with nonzero 
curvature. 
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1. Introduction 



Uchida, Sugitani and Ohya [H] studied a non-stratified airflow past a two-dimensional 
ridge in a uniform flow. Airflows around the ridge include an unsteady vortex shedding. 
Their study focused on airflow characteristics in a wake region. In general, before separat- 
ing from a boundary, the flow moves toward reverse direction near the boundary against 
the laminar flow direction (see [HI E] for example). There are several results related to 
the wake region in pure mathematics. Using the Oseen system is one of the mathematical 
approach to analyze the wake region. For the detailed discussion of the Oseen system, 
we refer the reader to [1]. In a convex obstacle case, the character of the system is el- 
liptic in front of the obstacle. To the contrary, its character changes into parabolic type 
(wake region) behind the obstacle (see [9] for example). Maekawa [12] considered the 
two-dimensional Navier-Stokes equations in a half plane under the no-slip boundary con- 
dition. He established a solution formula for the vorticity equations and got a sufficient 
condition on the initial data for the vorticity to blow up to the inviscid limit (see also 
[4]). His observation suggests that if the Reynolds number is high, the boundary layer 
immediately appears and the high vorticity creation occurs near the initial time and the 
boundary. Ma and Wang [11] provided a characterization of the boundary layer separa- 
tion of 2-D incompressible viscous fluids. They considered a separation equation linking 
a separation location and a time with the Reynolds number, the external forcing and the 
initial velocity field. In the Dirichlet boundary condition case, which corresponds to the 
boundary layer separation, the above mentioned work of Ma and Wang provides detailed 
information on the flow transition near the critical time. However, none of the above 
studies has shown the mechanism behind the reverse flow phenomena rigorously. In this 
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paper we observe a mathematical clue to analyze such reverse flow phenomena. More 
precisely, we consider the non-stationary Navier-Stokes equations with an initial datum 
having a parallel laminar flow (we define rigorously later). We show that the direction of 
the material differentiation is opposite to the initial flow direction and effect of the mate- 
rial differentiation (inducing the reverse flow) becomes bigger when the curvature of the 
boundary becomes bigger. Before showing such result, we show the parallel laminar flow 
cannot be a stationary Navier-Stokes flow near a portion of the boundary with nonzero 
curvature. 

Now we formulate our results in the following mathematical setting. Let Q be a either 
bounded or unbounded domain with smooth boundary in M 2 . The stationary Navier- 
Stokes equations are expressed as: 

(1.1) — vAil + (u ■ V)u = —Vp and div u = in £1 

with u\sn = 0. We need to handle a shape of the boundary dQ precisely, thus we set a 
parametrized smooth boundary tp : (— oo, oo) h-» IR 2 as |5 a <p(s)| = 1, 9(d s (p(s)) is a decreas- 
ing function, U_ 00 < s<00 (y9(s) C dQ, where 9{w) is defined by a vector w = r(cos sin #), 
9(w) := 9. For a technical sense, we need to assume that there are s and s 2 (we set Si = 0, 
so < $1 < S2) s.t. <p(sq) = (0,0), 9(d s <p(s))\ s =s = 0, \d%<p(s)\ is monotone increasing for 
s G [so, si], |<9 2 y?(s)| = 1/5 for s G [si, S2], where 1/5 is a constant curvature of a part of 
obstacle boundary U sg [ Sl) s 2 ]V?(s). Note that there are s and x, then <p(s) (s G [si, S2]) can 
be expressed as 



Here we mainly consider laminar type flows (we define the "laminar type flow" later) 
in a localized region near the boundary, thus we need to assume m / (no stagnation 
point) in Q near U§ 0<s<S2 </?(s). Boundary layers appear on the surface of bodies in viscous 
flow because the fluid seems to stick to the boundary dfl. Right at the surface the flow 
has zero relative speed and this fluid transfers momentum to adjacent layers through the 
action of viscosity. We need to express such phenomena in pure mathematics. To do so, 
we need to give the following coordinates. 

Definition 1.1. (Normal coordinate.) For s G [sq 5*2)5 let 

$(s, r) = $ v (s, r) := (9 s y?(s)) ± r + <p(s). 

We define _L as the upward direction. 

Remark 1.2. (An infinite channel with uniform width.) We give a typical example of 
Q. Assume that there is s 3 (> s 2 ) s.t. \d^ip(s)\ is monotone decreasing for s G [s 2 , s 3 ] and 
\8%(p(so)\ = |<9 2 <£>(s3)| = 0. Assume moreover that Q is a straight pipe except near the 
origin, namely, we assume that (p(s) = (s, 0) for s < 0, 9(d 8 cp(s)) is a constant for s > S3 
and —it < 9(d s ip(s3)) < 0. Then we find an infinite channel with uniform width as 



Q = {$( S) r) : s G (-00, 00), r G (0, 1)}. 
In this case, there are several results related to the Poiseuille flow (see [9 J for example). 



(1.2) 
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Remark 1.3. By a scaling argument with a cone shape (11.21) . we see that for any s±, s 2 € 
[si, s 2 ] with si < s 2 , 

I U% S1 $(^r0| = I U- =S1 $(,',r 2 )|. 

In order to state main theorems, we need to give several definitions. 

Definition 1.4. (Normalized streamline and its direction.) Let 7 a be in Q near U So<s <j 2 (^(s) 
which satisfies 

<9 s 7a(s) = \ T~n (7a0)), 7a(0) = a e near Us 0<s<52 ip{s). 

Moreover, we assume that the flow is to the rightward direction (laminar flow direction), 
namely, 

0<(d 8< p(s),^(s,r))\<l 

for s e [so, s 2 ] and sufficiently small r, where (•, •) means inner product. 

We have already defined that there is no stagnation point near the boundary. Thus the 
above definition is well defined. 

Definition 1.5. (Poincare map.) For fixed s and Si sufficiently close to each other, let 
Smin be the minimum of s' > for which there exists r = r(s') such that $(si,r(s')) = 
7<&0,r)(s')- Let L(r) = L s>si (r) = r(s min ). 

Definition 1.6. (Classification of laminar type flows.) The velocity u near the boundary 

U s ' e [s 0i s 2 ] l f{s') is called: 

"Strong diverging laminar flow" iff L(r)/r > C > 1 near the boundary, 
"Weak diverging laminar flow" iff L(r)jr > 1 and R(r)/r — > 1 (r — > 0), 
"Parallel laminar flow" iff L(r)/r = 1 near the boundary. 

Remark 1.7. In the parallel laminar flow case, we see that |w(7 a (s))| = \u(a)\ if a and 
7a (s) are sufficiently close to U s€ [s 0t s 2 ]f(s). 

In this paper we only consider the parallel laminar flow case. In the parallel laminar 
flow case, we see that 

(1.3) 9{d sl{0 . x2)i s ))\s=o — for sufficiently small 22 > 0. 
In this case, it is natural to set 

(1.4) 1*1(0, x 2 ) = a.\Xi — ^-xl(=: h(x 2 )) and w 2 (0,a; 2 ) = («i,a 2 >0) 

in order to satisfy (II. ip for sufficiently small x 2 > 0. In this case, we can say that the 
boundary layer thickness is cti/a 2 . This setting is based on "Poiseuille flow" (see Remark 
11.21 and [8 J for example). Thus the parallel flow satisfying ( II .4p is one of the candidate of 
the solution to the stationary Navier-Stokes equation. 

Remark 1.8. We see from Remark 11.71 that 

(1.5) h(x 2 ) = |ix($( 5 ,x 2 ))| 
if $(s,x 2 ) is sufficiently close to U se [ S0j s 2 ]y?(.s). 

The first main result is as follows: 
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Theorem 1.9. (Stationary case.) Let us choose v > arbitrarily and fix it. Then there 
is no parallel laminar flow satisfying the stationary Navier-Stokes equation U.l\) and ( TTjU 
near U%^ip(s). 

Next we give a result in the non- stationary Navier-Stokes case. The non-stationary 
Navier-Stokes equations are expressed as 



;i.6) 



d t u — uAu + (u ■ V)u = — Vp, div u = in QcR 2 
u\dn = and u\ t=0 = u . 



Throughout the non-stationary case we assume that p(x, t) \ t= o is smooth in space valuable. 
If Uq is smooth and in L 2 (Q) (which means finite energy), then the solution (u,p) is a 



global- in-time smooth solution (see [10]). By Theorem ll.9[ we see that if we take an 
initial datum which has a parallel laminar flow near the boundary, the solution u must be 
a non-stationary flow in some small time interval near the initial time. In this case, we 
observe a mathematical clue to the separation phenomena. More precisely, we see that 
the direction of the material differentiation is opposite to the initial flow direction and the 
effect of the material differentiation becomes bigger when the curvature of the boundary 
becomes bigger. Let be more precise. The initial datum having a parallel laminar flow 
satisfies 

(1.7) (^^(x),d s (p(s)\ 1 as x ->■ ip(s) (s G [si, s 2 ])- 

At some time t > 0, if the solution u(x,t) has the reverse direction near the boundary 
against the laminar flow direction, we should have 

(1.8) ^- — -(x, t), <9 s <^(s)^ — > — 1 as x — > (p(s) (s G [si, S2]) for some t > 0. 

To obtain such changing direction property, we can expect that the solution should satisfy 

(A M (M)^(.)) ^ 
{uo(x),d s (p{s)) 

for a long time interval up to the time t, where D t u := dtu + (u ■ V)« is the material 
differentiation. We show that the solution u(x,t) to the non- stationary Navier-Stokes 
equation satisfies (11 .9p at the initial time, if the initial datum has a parallel laminar flow 
near the boundary. The following is the second main theorem. 



Theorem 1.10. Assume that Vp(t,x)\t=o is smooth. If the initial datum satisfies \l-4\) 
and has "Parallel laminar flow" , then we have 

(A^, t )U.^M> = _^_. <0 

x-+tp(s) {u {x),d s ip(s)} oai b l 

Recall that the value ol\Ioli is defined as the boundary layer thickness. In general, the 
pressure is a non-local operator. However the values in the above theorem are depending 
only on the behavior of the flow near the boundary U sE \s 1} s 2 ](p(s). 

Remark 1.11. There are direct and indirect evidences for the validity of the "Kutta 
condition" in restricted regions (see [2]). The method used in the above theorem may 
give another support for the validity of the Kutta condition in pure mathematical sense. 
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Remark 1.12. For the Poiseuille flow case, the pressure is depending on v. More pre- 
cisely, if 

fi = {iGl 2 :fl<i 2 <l} 1 

then u = [xi — x 2 ,0), p = (u/2)x2 is one of the stationary solution to (11.11) . To the 
contrary, for the whole space case (without any boundary effect), the pressure (in the 
sense of "mild solution" ) is expressed as 



P = R i R 



where Rj is the Riesz transform (see [5j [7J H3]). In this observation, the dependence on 
the viscosity appears due to the Dirichlet boundary. 



2. Proofs of the main theorems 

Proof of the first main theorem. First we estimate Au. For s G [si, S2] that is selected, 
we call the point <f(s) to be Q. That is Q = <p(s). Then, we use the orthonormal frame 
e 1 = d s ip\s, e 2 = (dsip)- 1 -^ at the point Q to construct a cartesian coordinate system with 
the new Xi-axis to be the straight line which passes through Q and parallel to the vector 
ei, and the new X2-axis to be the straight line which passes through Q and is parallel to 
the vector e^- Then, for the given vector field u in an open neighborhood near dQ, we 
define for each y 6 fi near dfl, the two components of u(y) with respect to the e\ direction 
and the e 2 direction. It is expressed as 

u(y) = v 1 (y)e 1 + v 2 (y)e 2 , 

where v 1 (y) := (u(y), e\) and v 2 (y) := (u(y), e2). Since <p(\si, s 2 ]) is known to be a circular 
arc with radius 5, let C to be the point at which the center of the circular arc (^([si,^]) 
is located. Let (s, r) to be the coordinate representation of the point y in the coordinate 
system based at Q which is specified by the orthonormal frame {ei, €2}- That is, the point 
y can be realized as y = Q + se\ + re 2 . Since the vector field u is a parallel laminar flow 
near dQ, the streamline of u which passes through y should be a circular arc with radius 
\Cy\ and centered at C also. Here, \Cy\ is length of the line segment Cy That is, the 
distance between C and y. Moreover, we have 

\u(y)\ = h(\Cy~\ - 5) 

with 6 to be the radius of the circular arc £2]), and h{r) = a±r — ^-r 2 . So, in order 

to compute |w(y)|, it is enough to compute \Cy\ by means of geometry. Let L x to be the 
straight line which passes through Q and is parallel to the e\ direction. That is L x is the 
new x — axis of the new coordinate system. First, let A to be the unique point on L x such 
that the line segment QA is perpendicular to Ay. Let D to be the point of intersection 
between the line segment Cy and the line L x . Observe that \Ay\ = r, \QA\ = s, and 
\CQ\ = 5. We first compute \QD\ and \DA\ through the following observation. Since 
ADAy and ADQC are similar triangles, it follows that 

\DA\ = \QD\ 
r 5 
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Hence, by substituting \DA\ = (r/5)\QD\ into the equation s = \QD\ + \DA\, it follows 
that 

T til 

\DA\ = -— and \QD\ 



5 + r 5 + r 

This further gives 



2 12 f / \ 2 1 2 



Hence, we have 

\C^\ = m + \DC\ = {(5 + r) 2 + s 2 }l 
As a result, we have \u(p) \ = h{{{5 + r) 2 + s 2 }^ — 5), and that 

v 1 = h({(6 + r) 2 + s 2 } 1 ? - 6) ■ cos(%)) 
v 2 = h({(5 + r) 2 + s 2 }5 - 5) ■ sin(%)), 

where, 6(y) is defined by the inscribed angle between Dy and Ay. Since ADAy is a right 
angled triangle, it follows that 



sin(%)) = -sin(|%)|) = -^ 



(210) \Dy\ {(s + ry + s^ 

cos(%)) = co S (\9(y)\) = — * + r 

{(o + r) 2 + s 2 \2 

Hence, it follows that 

/ = v\Q + se 1 + re 2 ) = h{{{5 + r) 2 + s 2 }^ - S) ■ * ' 



(2.11) + + 

v 2 = V 2 (Q + S ei + re 2 ) = -h({(6 + r) 2 + s 2 }^ - 5) 



i • 
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{(5 + r) 2 + s 2 } 
Therefore we obtain, through direct computation that, 

dlv 1 = —a 2 , dlv 2 = 0, d 2 v 1 L =0 = — ^— : \ol\ — a 2 r — ^- ) and d 2 v 2 = 0. 

r + o \ r + o J 

From (11.51) . we see «(7$(o,r)( s )) = r)) = /i(r). Thus the material differentiation can 

be calculated as 

((« ■ V)u) (7*(0,r)(s)) = <9 S ( M (7<£(0,r) (■§))) | s =g 

u 

= |M(7*(0,r)(s))|9 s (|-|(7*(o,r)(s)))| s = 5 



+<9 S K7*(0,r)(s)) | ( |-r(7*(0,r) (s))) | 



|w(7*(0,r)(s))|9«(y (7*(0,r)(s)))|* 

--^(r)(«9^( S -). 



Note that the above decomposition (into "curvature part" and "acceleration part") is 
already done in PQ. Therefore 

(2.12) i/Att($(s,r)) - ((it ■ V)u) ($(s,r)) = P(r)^(s) + P ± (r)(9 s ^) ± (s) 
for s G [si, S2], where 

(2.13) P ( r) = Jj!^-^- T W-a 2 ), P\r) - *< r > 



r + 5 r + 5 (r + 5) 2 /' 1 r + 5 

In order to see the pressure precisely, we need the following "normalized pressure-line" 
and "boundary of the level set of the pressure" . Let 

d s q a (s) = |^| (?-(*)). (|^|) (i"(r)) for |V P |^0 

with g a (0) = a and ^(0) = a. Moreover let us define a re-parametrized q as follows: 

Definition 2.1. (Poincare map on the pressure lines.) For sufficiently small e > 0, let 
Smin be the minimum of s' > for which there exists r = r(s') such that 

?#( 8 ,r)(s') = qi(s+e,r)( T ( S '))- 

Let 77(e) = r)(e,s,r) = q<z>( s>r )(s min ). We denote U^(e') as U e e , =0 r)(e' , s, r). 

Now we consider the level set of p precisely. We choose the following finite points: 
{r7(e'(0)),r/(e'(l)),r/(e'(2)),--- ,T](e'(N - l)),r](e'(N))} in order to satisfy that e'(0) = 0, 
e'(iV) = e and \r)(e?(£)) — rj(e f (£ — 1))| (£ — 1, • • • N) are the same value. Then we set Vp 
(in the discrete setting) from the level set of p 

P[l{ [))) ~ \v(e'(£))-v(e'(i + l))\ \v(eW-v(e'(£ + l))\ 
and let t be a tangent vector defined as 

n , m) = V(e\£)) - r,(e\£ + I)) 

1 1 )] W(£))-v(e'(£+i))\' 

Then we have in the discrete setting (piecewise linear approximation): 

N 

/ Vp-f = lim J2Vp(e'(£))-t(6'(£))\r)(6'(£))-r)(e'(£ + l))\ 
= p(ij(e,s,r)) -p($(s,r)). 

Therefore 
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(2.14) |Vp($(s,r))| = lim 



Vp-t 

U e Jj( £ ') 



I U e 77(e')| 

Since the direction of Vp, namely is determined by and jf^^j , and u has a 
parallel laminar flow, we see that 

/ JVp_ d s $(s,r) 

\\Vp\ [q>[S > r)) '\dMs,r) 
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is independent of s. Thus there are positive s(r) and f(r) (independent of s) such that 



p(<t>(s, r)) = p(?i ( . >r) (0)) = p(gi (SlP) (r(r))) = p(S(a - s(r), 0)) 



(a,r)l 



and 



Since 



p(77(e, s, r)) = p(g^ (a+eir) (0)) = p(g* (a+e , r) (r(r))) = p($(s + e - s(r), 0)). 



Vp($(s, 0)) = i/Au($(s, 0)) = - a 2 )0 s¥ > 

o 

derived from the boundary condition, we see 



|p($(s + e - s(r), 0)) - p($(s - s(r), 0))| 



Vp($(s',0)) -<p(s') ds' 



o 



ei/ — - a 2 \ 
o 



Thus we have the following pressure estimate: 



|Vp($(s,r))| = hm ' ' , 



e^O | U £ rj(e' 



On the other hand, since 9 ((d s $)(s + e)) > (<9 s g<j>( s ,r)(s + e)) for any e, there is a 
negative number r = r(e, s,r) < such that <2$( s+e r ) (t) = r)(e,s,r). Thus -D e := 

V') is well defined and \im e ^. Q D e /\D e \ is a right triangle since 



l \ e l l u 



-T(e',s,r) "*(s+e',r) 



((<9 s $)(s, r), ((9 e ?7)(0, s, r))) does not change for any e > and 

(^9*(«+er)( T )' (deV)( e i s > r )) = f° r an y e > 0. It means that angles of the three corners 
do not change for any e and one of them is tt/2. More precisely, among the three sides 
of the triangular region D e , the side Uo< e '< e r](e' , s, r), which is a part of the pressure line 
passing through $(s,r), is perpendicular to the level set of the function p, which passes 
through $(s + e,r), a£ the point of intersection r)(e,s,r). In other words, we have the 
following observation. 

• the right angle | of the triangle D e is located at the vertex r](e,s,r), which is 
the point of intersection between the side Uo< e '< e ^(e / , s, r) and the level set of p 
passing through $(s + e, r). 

So, the arc segment | U s < s '< s+e $(s',r)| should be the longest side of the triangle D e 
whose opposite angle in D t is the right angle | located at the vertex r)(e, s, r) of D e . 
According to the above reasoning, it is not hard to see that 



(2.15) 



\P(r) 



U < € '<eV(e',s,r) 
[(P(r)) 2 + (P^(r)) 2 ]5 ~ eTo+ | U s < s ,< s+e $(s',r) 



= lim 



In accordance with the relation in (12.151) . it follows that 



e e | U s<s , <s+e $(s',r) 

hm ; — = hm 



£-^0+ | U < e '< £ T]( e ', S , r ) I I ^s<s'<s+e ®(S', r) I I U < e '< e T]( e ' , S i T ) I 

e \(P(r)) 2 + (P^(r)) 2 ]l 

= (^\U s < s ,< s+e <S>(s>,0)\) ' \P{r)\ 

[(P(r)) 2 + (P ± (r)) 2 ] 1 2 



\P(r)\ 

5 \ [(P(r)) 2 +(P ± (r)) 2 ]^ 



\5 + rj \P( r )\ 

Hence, it follows from relation ( 12. 14ft and the above that 

ivy i»i l ( M [(P(r)) 2 + (P ± (r)) 2 ]^ 

|Vp($(s,r))| = - ar 2 | ' 1 



5 'yS + rJ \P(r)\ 

However, since we naturally have |Vp($(s,r))| = [(P(r)) 2 + (P _L (r)) 2 ]^, it follows from 
the above relation that 



5 \5 + r J \P(r) t 

nr-oocmn P ( >r\ — i,i _ 

■ r+S (r+S) (r+5) 2 



The above expression, together with the expression P(r) = i^(-^h — t^ttt — , h j_ r ) )2 — oc 2 ) 



directly implies that 
(2.16) \—-a 2 \ 



5 \5 + r 



«i a 2 r h(r) 

a 2 



r + 5 (r + 5) (r + 8)' 



Our task here is to derive a contradiction from the relation in (I2.16p . First, let us deal 
with the case in which I ^- — «2 1 > 0. In this case, since (;qjg — ~ ( r +g)3 ~ a 2) — Q2, 

as r — >■ + , it follows that — r*+jn — — a 2 ) will have the same sign as that of 

— «2, for all sufficiently small r > 0. That is, for all r > to be sufficiently small, the 
following assertion holds. 

(7+6 ~~ 17+5) ~~ (tSf ~~ " 2 ) is P ositive ^ and omv if x ~ " 2 is P ositive - 
As a result, the relation in (12 . 1 6[) can be rephrased as 

( 2 - 17 ) T ~ ft 2 ^— = — -7 ~ 7—^ - 7- , ^ - a 2 . 



5 J \5 + rJ r + 5 (r + 5) (r + 5) 2 

However, by some algebraic computation, relation as given in ( I2.17P is equivalent to saying 
that the following relation holds for all r > sufficiently close to 0, 

5_ \ h(r) a 2 r 



5 + r J (5 + r) 2 r + 5' 

However, the above relation says that, for r > to be sufficiently small, the left hand 

side oi 2 [j^p — 1], which is strictly negative, equals the right hand side + 7^§y, which is 

strictly positive. This is absurd, and hence we have arrive at a contradiction in this case. 
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The case of ^ — a 2 = is even easier to handle, since in this case, relation (I2.16P 
immediately gives 

a.\ a 2 r h(r) Sa 2 h(r) 

~5 + r r + 5 (5 + r) 2 " 2 ~~ 5 + r (5 + r) 2 " 2 

which eventually leads to the relation — 1] = J^) 2 + which is again absurd. 

So, we also have a contradiction. 

Proof of the second main theorem (Non-stationary case.) Let p(x) := p(x,t)\ t=0 and 
u(x) := u(x,t)\ t=0 . From (I2.13p . we see that 

Thus we have 

(2.18) Vp($(s,0)lr=o = z/Aw($(s,r))| r=0 = f (y - <* 2 ) 

derived from the Dirichlet boundary condition. In order to estimate Vp only using smooth- 
ness and (I2.18p . we use "normalized pressure- line" and "boundary of the level set of the 
pressure" defined in the stationary case. The key is to estimate | U^/ =0 C( e ')l ( we define it 
later) which is similar to | U e e , =0 rj(e')\ defined in the stationary case. 

Definition 2.2. Here we give three definitions which are needed in a geometric obser- 
vation. 

• Let s = s(s, r) and f = f(s, r) be such that 

<li(s)(r) = <5>{s,r) for s G [si,s 2 ]. 

• For sufficiently small e > , let s m i n be the minimum of s' > for which there 
exists r = r(s') such that qF(s') = Q^ s+e )( T (s')), where F = Q^f s \(r(s,r)). Let 
C(e) = C(e, s,r) = q F {smin)- We denote U e C(e') as U^, =0 C(e', s, r). 

• Let s = s(e, s, r) and f = f(e, s, r) be such that £(e, s, r) = f). 

Due to the smoothness of p and (I2.18p . we can estimate s and s as 

(2.19) s — cr 2 < s < s + cr 2 and s + e — cif 2 < s < s + e + C2r, 

where c is a positive constant independent of e, s, r and f. Let us set D tr (which is 
similar to D e in the stationary case but different) as follows: 

A,r :=U«, =0 Up =0 C(e',s,r'). 
Since all of the four angles are tt/2, 

lim — 

e,7 — >0 |£) I 

r/e=const. ' e ' r ' 

is a rectangular. Thus we can see that for sufficiently small e > 0, there are R, R and i? 
such that 

(2.20) (1 - e)r < r < (1 + e)r 
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for r < R with r/e = const., r < R and r < R. From (I2.19p . we have 



and then 



s + e - c(l + e) V <s<s + e + c(l + e) V 



e — cr 2 — c(l + e)r 2 < s — s < e + cr 2 + c(l + e)r 2 . 



Due to the smoothness of p and (12 . 18[) , we also see that 



1 - cr 2 < 



( VP 

Vivp| 



($(s,r)),9^(s)> < 1. 



By using a piecewise linear approximation, we can estimate | U e C( e ')l as 



N 



U e C(e')l 



lim / 



As, 



where Ar = 0, Ar fc = tan 6»($(s + (k - l)As,r + 1 Ar fc /)) for fc > 1, 6»($(s,r)) is 
defined as the angle between (Vp/| Vp|)($(s, r)) and d s ip(s), and we choose As in order 
to satisfy 



iVAs 



5 



s . 



Note that As tends to as N — > oo. It means that Ar^ tends to as iV -)• oo for fixed 
fc. Then we see that 



lim \ Arfc + 

n— >oo,iV— >-oo ' 
fc=0 



r. 



By the above estimates, we have the following lower and upper bound of | U e C( e ')l : 



(2.21) 



U e C( e ')l > (l-e)^(e-cr 2 -c(l + e)V), 



U e C(6')| < 



1 + e r + 5 



(e + cr 2 + c(l + e) 2 r 2 ) 



1 — cr 2 5 

for sufficiently small e and r with e/r = const. Note that we can take e > arbitrarily 
small. We have from the explicit representation of the gradient of the pressure (here we 
consider the case aci/5 — a 2 > 0): 



(Vp(<p(s)),d a <p(s)) = \Vp(<p(s))\ 

= |Vp($(s,r))| 
1 



Vp 



(<p(s)),<p(s) 



lim 



lim 



o|U e C(e')l 
1 



Vp(x) ■ t dx 



— \p(<p(s + e)) -p(<p(s))\ 



o|U^C(e')l 

,. eu(ai/5-a 2 ) 
— hm — ; — — 7-: — . 

e ->o |U e C(e')l 

Even in the case a±/5 — a 2 < 0, we have the same estimate as above. Recall that 

an a^r /i(r) 



(uAu, dip(s)} = v 



r + 5 (r + 5) (r + 5) 2 
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«2 



and 



(u ,d s ip(s)} = a x r - yr 2 . 



We already have the lower and upper bounds of | U e C( e ')l ( see (12.211) ). we have that 
{D t u, d s (f(s)} ((uAu - Vp),d 3 (p(s)) 
(u ,d s ip(s)) (u ,d s ip(s)) 

Um (i/Ati($(a,r)),0,p( a )) - (eKgi/j - a 2 )/| U e CfeQD 

(u ($(s,r)),<9 s y?(s)) 

e/r=const. 

voti v 



5cti 5 2 

for s G [si, S2]. Therefore we have the the desired estimate. 
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